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The negative pressure accompanying gravitationally-induced particle creation can lead to a cold 
dark matter (CDM) dominated, accelerating Universe (Lima et al. 1996 [lj) without requiring the 
presence of dark energy or a cosmological constant. In a recent study Lima et al. 2008 @] (LSS) 
demonstrated that particle creation driven cosmological models are capable of accounting for the 
SNIa observations of the recent transition from a decelerating to an accelerating Universe. Here 
we test the evolution of such models at high redshift using the constraint on z eq , the redshift of the 
epoch of matter - radiation equality, provided by the WMAP constraints on the early Integrated 
Sachs- Wolfe effect (ISW) Q. Since the contribution of baryons and radiation was ignored in the 
work of LSS, we include them in our study of this class of models. The parameters of these more 
realistic models with continuous creation of CDM is tested and constrained at widely-separated 
epochs (z ~ z eq and z ~ 0) in the evolution of the Universe. This comparison reveals a tension 
between the high redshift CMB constraint on z eq and that which follows from the low redshift SNIa 
data, challenging the viability of this class of models. 

PACS numbers: 98.80.-k, 95.35.+d,95.30.Tg 



I. INTRODUCTION 

A large amount of complementary cosmological data have established the surprising result that the Universe expe- 
rienced a recent transition from a decelerating to an accelerating expansion [5J] . The current data are consistent with 
a flat, Friedman-Lemaitre cosmology whose dominant components at present consist of matter (cold dark matter plus 
baryons) and a cosmological constant (A) or, equivalently, the energy density of the vacuum; the ACDM concordance 
model [6]. The matter and radiation dominate earlier in the evolution of the Universe, leading to its decelerating 
expansion, while the cosmological constant or vacuum energy component has come to dominate more recently, driving 
the accelerating expansion. While the concordance model provides the simplest, most economical explanation which 
is consistent with all extant data, provided that it is fine-tuned to fit that data, it is by no means the only candi- 
date proposed to explain the accelerating expansion of the current Universe. Dozens (if not hundreds!) of alternate 
possibilities have been explored in the literature (see Lima et al. 2008 ,2j (LSS) for an extensive list of references) 
and many of them remain viable candidate models. Since the astronomical community is planning a variety of large 
observational projects intended to test and constrain the standard ACDM concordance model, as well as many of the 
proposed alternative models, it is timely and important to identify and explore a variety of physical mechanisms (or 
substances) which could also be responsible for the late-time acceleration of the Universe. 

In principle, any homogeneously distributed "exotic" fluid endowed with a sufficiently negative pressure, p < — p/3, 
so-called dark energy, will suffice. The simplest candidates for the dark energy are a positive cosmological constant 
A, or a non-zero vacuum energy density. In both cases the equation of state of the exotic fluid is w = p/p = — 1. 
All current observational data are in good agreement with the cosmic concordance (ACDM) model, consisting of a 
cosmological constant or vacuum energy (w\ = — 1) whose current density parameter is Qa w 0.74, plus cold dark 
matter and baryons (u>m = 0) with 17m = ^cdm + = 1 — ~ 0.26. Nevertheless, the ACDM model encounters 
several challenges. For example, the value of the energy density of the vacuum required by the data (f^A ~ 0.26) 
needs to be fine-tuned, exceeding the estimates from quantum field theories by some 50 - 120 orders of magnitude. If 
the vacuum energy (or A) differs from zero, why does it have the value needed for consistency with the data? 

LSS noted recently 0] that the key ingredient required to accelerate the expansion of the Universe is a sufficiently 
negative pressure, which occurs naturally as a consequence of cosmological particle creation driven by the gravitational 
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field [7|. While completely different from it, there are some analogies between this class of models and the continuously 
accelerating, never decelerating, expansion of the steady-state cosmology driven, in the Hoyle-Narlikar Q version, by 
C-field induced particle creation. The late-time evolution of the class of models explored by LSS is qualitatively very 
different from that of the standard, ACDM concordance model in that they lack a cosmological constant and, when 
baryons are included, the ratio of the dark matter to baryon densities evolves with time or redshift. Nonetheless, the 
LSS [2| analysis of the late-time evolution of such a CDM-dominated model with no cosmological constant or dark 
energy, lacking baryons and radiation, showed that it is quantitatively possible to account for the SNIa data which 
reveal a recent transition from decelerated to accelerated expansion. Indeed, while the flat ACDM concordance model 
provides a good fit to the SNIa data, the best fit actually corresponds to a positively curved (closed) ACDM model, 
so that for some choice of parameters the models considered by LSS may provide an even better fit to the supernovae 
data. For a subset of the models explored by LSS to be discussed here, the early-time, high redshift (e.g., z > 10) 
evolution is indistinguishable from that of the ACDM concordance model, so that the parameters controlling the 
early evolution of these models must agree with those of the concordance model. To set the stage for the constraints 
and tests discussed here, we note that for the ACDM concordance model with = 1 — Qm = 0.76, the redshift of 
the transition from decelerated to accelerated expansion is Zt — 0.79, the redshift of the epoch of equal matter and 
radiation densities (related to the early, Integrated Sachs- Wolfe (ISW) effect) is z eq — 3223 and, the present age of 
the Universe in units of the Hubble age (H^ 1 ) is Hota = 1.00. 

Since at present baryons are subdominant to CDM and the contribution of radiation is negligible, the neglect of 
baryons and radiation by LSS was not a bad approximation for their study of the recent evolution of the Universe. 
Here, we revisit the LSS model of dark matter creation driven acceleration in a more realistic model incorporating 
baryons and radiation. While it will be seen that there is a subset of models with baryons and radiation whose late time 
evolution (i.e., H = H(z) for z < 5) is sufficiently close to those studied by LSS and to the ACDM concordance model 
(and, therefore, they too will be consistent with the SNIa data), these models differ from those studied by LSS in their 
early evolution, in particular at the epochs of recombination and of equal matter and radiation densities. To further 
test these models we impose an additional, high-redshift constraint derived from the early Integrated Sachs- Wolfe 
(ISW) effect, provided by the WMAP Q value of the redshift of matter - radiation equality, z eq = 3141 ± 157. 

In § [II] we describe the class of models under consideration and solve for the evolution of the Hubble parameter 
as a function of redshift (H = H(z)). In § III we compare the predictions of this class of models as a function of 
the two parameters related to the particle creation rate to the SNIa data [|| and to the constraint on the redshift 
of the epoch of equal matter and radiation densities provided by the WMAP observations of the CMB temperature 
anisotropy spectrum 0] ■ Our results are summarized and our conclusions presented in § IIV1 



II. PARTICLE CREATION DRIVEN COSMOLOGY 



Consider a sufficently large comoving volume (U), representative on average of the Universe which, at some time, 
contains N = nV CDM particles (wcdm = 0) . If CDM particles are being created at the expense of the gravitational 
field 0, then 

^ = r, (i) 

N dt K ' 

where T is the creation rate (number of particles per unit time), assumed to be uniform throughout the Universe. In 
this case, since the mass density in the CDM is proportional to the number density, 

d[\n(p CBM V)} = r 
dt 



so that 



Pcdm = Pcdm,o(1 + z) 3 exp[- / Tdt']. (3) 



In eq. 3, z is the redshift (V = Vb(l + z) 3 ) and the subscript is for quantities evaluated at the present epoch 

The models we consider are defined by the choice for the particle creation rate, T. The most natural choice would 
be a particle creation rate which favors no epoch in the evolution of the Universe, such as T = 3/3 H, where H is the 
Hubble parameter and the free parameter (3 is positive. However, in the absence of baryons and radiation, it is easy 
to show (see LSS) that such a Universe will always decelerate if (3 < 1/3 and will always accelerate if (3 > 1/3. The 
inclusion of baryons and radiation, whose effect is to drive a decelerated expansion, opens the possibility of a transition 
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from an early, matter-dominated deceleration, to a late, CDM-driven accleration. In this case the transition from 
decelerating to accelerating expansion is set by the baryon density and it occurs "naturally" , late in the evolution of 
the Universe at low redshift, without any need for fine-tuning. However, it will be seen that the early evolution of 
models with this form for the particle creation rate (T oc H) are inconsistent with the early ISW data. This problem 
may be alleviated by adding a second creation term for the particle creation rate which is "tuned" to the current 
epoch, r oc Hq 7 where Ho is the present (to) value of the Hubble parameter. As a result, the class of particle-creation 
driven models to be studied here are defined by, 

r = 3jH + 30H, (4) 

where < {7, 0} < 1. In this case, the CDM mass/energy density evolves as, 

Pcdm = Pcdm,o(1 + 2) 3(1_/3) exp[37(r - to)], (5) 

where r = Hot and To = -Ho^o are the age and the present age of the Universe in units of the Hubble age (Hq 1 ). For 
Ho = lOO/i kms^Mpc" 1 , iJ _1 = 9.78k- 1 Gyr. In our analysis the HST Key project result @ h = 0.72 ± 0.08 is 
adopted, so that Hq 1 = 13.6 ± 1.5 Gyr. 

For simplicity, as is the case for the cosmic concordance model, our considerations are limited to flat cosmologies, 
so that Qm — f^B + ficDM = 1, neglecting the very small contribution from the radiation density at present, when f2e 
and f^cDM are evaluated. For flat models, the general Friedman equation reads, 

="R( 1 + z ) 4 + r! B(l + ^) 3 + (l-r!B)(l + ^) 3(1 - /3) exp[37(T-To)]. (6) 

In comparing the model predictions with the SNIa data, the radiation density term (SIr < 10~ 4 ) may be safely 
neglected. 

At the time of Big Bang Nucleosynthesis (BBN) the Universe is radiation dominated and only the baryon density 
(along with the radiation density) plays an important role. BBN cares about neither the dark matter or the cosmo- 
logical constant for the ACDM model, nor about the creation of dark matter for the {/3, 7} model considered here. 
However, since the creation of dark matter is accompanied by a "creation pressure" p c = —Tpdm/3H (see LSS), and 
the evolution of the dark matter density differs from the "usual" (1 + z) 3 evolution for conserved particles, the late 
time growth of perturbations in this model will likely depart from that in the concordance model. At very early times 
near the epochs of equal matter and radiation densities and recombination, w c = p c /pr>M —> —ft- So if, as will be 
seen below, (3 = or f3 <C 1 are favored, the early growth of perturbations in the {/3, 7} model will track that of the 
concordance model. As a result, for the second observational constraint we require that the model-predicted redshift 
of the epoch of equal matter and radiation densities, 1 + z eq = pm/pr, agree with that determined by the WMAP 
observations of the early ISW effect Since the ACDM model is consistent with these observations, the early-time 
(i <SC £0) evolution of these {/3, 7} models must track closely that of the standard ACDM model. 

Before considering the general case, where {fte, P, 7} 7^ 0, it is instructive to examine several simplified cases where 
one or more of these parameters is set equal to zero. 



A. The LSS Model: Q B = 



Ignoring the contributions from baryons and radiation, this is the model studied by LSS 2]. In this case the 
Friedman equation simplifies to, 



H 



1 / da\ 
a \d,T I 



i 2 



= (l + z) 3 ( 1 -«exp[3 7 (T-To)], 
which has a solution for the scale factor a(t) (or the redshift, z) as a function of time, 



1 



£(*) = 

a (I + 2) 



1-7-/3 



(e 3 ^/ 2 - 1) 



3(1-/3) 



Evaluating this at z = 0, where t = To and a = do, relates the present age, To, to the {/3,7} parameters, 

1-/3 

exp(3 7 T /2) = 



(7) 



(8) 



(9) 
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The present age of the Universe (t = to, r = tq) is 

tq = H t 

Solving for the time - redshift relation 

exp(37r/2) = 1 + 



2 -ln(- ^ 



37 '1-7-/3' 



1-7-/?, 

Using the z versus r and To relations above, the Friedman equation for H — H(z), reduces to the simple form 

H/H = (1 - /3)- x [7 + (1 - 7 - + zf^' 2 }. 



(10) 



(11) 



(12) 



It is easy to confirm (see LSS for the details) that if 7 = 0, the expansion of the Universe always decelerates for 
< (3 < 1/3 and always accelerates for 1/3 < (3 < 1. In this case there is no transition from an early decelerating 
to a late accelerating Universe. However, with 7 / (0 < 7 < 1) the redshift, Zt, of the transition from early-time 
deceleration to late-time acceleration is given by 



l + z t = [ 



27 



(l-3/?)(l-7-/?) J 



3(1-/3) , 



(13) 
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FIG. 1: The relative distance modulus (A(m — M)) - redshift relations for a variety of models with and without particle 
creation. For models without particle creation the ACDM model (solid, black) and the Einstein-deSitter model (solid, orange) 
are shown. Several particle creation models are shown, with and without baryons, for a selection of {/3, 7} choices. The inset 
figure shows the same curves along with the central values from the SNIa data 



As shown in LSS, there is a range of {7, (3} values whose fit to the SNIa data is at least as good as the fit provided 
by the parameters of the ACDM concordance model. In Figure [TJ we illustrate this by comparing the relative distance 
modulus - redshift relations for several combinations of {7, /?} with that for the ACDM model and, for comparison 
with a model which leads to a poor fit, the Einstein - deSitter model. Note that for z S 1-5, many of the {f3,j} 
models, with and without baryons, are indistinguishable from the ACDM model. While the curves for the best fit 
values for (3 = and {fie, 7} = {0, 0.63} and {0.042, 0.66} are very close to each other, the curve for 7 = 0, fie = 
0.042 and (3 = 0.61 deviates noticeably from them, especially at z ^ 1. Furthermore, the curves for (3 = and 
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FIG. 2: The 68% (dark blue) and 95% (light blue) contours in the j3 — 7 plane derived for J7b = from the SNIa data [jj, along 
with the 95% confidence band (black; see the inset which zooms in on this band) consistent with the early ISW constraint that 
Zeq = 3141 ± 157 4\. The SNIa data are best fit at /3 = and 7 = 0.63. The solid black curve, shown in more detail in the 
inset, is for the best fit relation between (3 and 7 = 7° (see eqs. 16,17) consistent with the observed value of z eq [J]. 
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FIG. 3: The likelihood functions for 7 = 7° inferred from the SNIa data (solid/gray) and from the z eq constraint (dashed/blue) 
for (3 = and Q B = 0. 



7 = 0.52 (0.50), corresponding to the z eq constraint with (without) baryons, deviates significantly from those curves 
which provide the best fits to the SNIa data. 

For the Union 2008 set of SNIa data Q adopted here, the best fitting combination of parameters (for Qb = 0.042) 
occurs for j3 — and 7 = 0.63, corresponding to zt = 1.26. While this value for the transition redshift may seem 
surprisingly high compared to z t = 0.79 for the ACDM concordance model with Qjvi = 1 — = 0.26, as may be 
seen from Figure Q] the relative distance modulus - redshift relations for these two models are very similar at redshifts 
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z ^ 1.5. For [3 = and 7 = 0.63, the present age of the Universe is tq = Hota = 1.05, corresponding to to — 14.3 Gyr. 

At high redshifts, including the contribution from radiation (but not yet from baryons), the Friedman equation is 
modified to 

so that the redshift of equal matter (CDM) and radiation densities in this class of models, without baryons (Ob = 0), 
is given by 

W + z eq f^={^-^ . (15) 

Fixing z eq , along with the observed value of Or, provides a complementary constraint on the 7 — (3 relation to that 
from the SNIa data. For Ob = 0, the 7(^3 = 0) = 7 — (3 relation is 

7 = (1 - - [0 R (1 + z e9 )] 1/2 (l + z eq fP' 2 }. (16) 

Accounting for the still relativistic neutrinos and for h = 0.72 ± 0.08, R = 8.065(1 ± 0.22) x 10~ 5 . For the 
combination of parameters which best fits the SNIa data [3j, = and 7 = 0.63, z eq = 1661, a redshift which is 
much to small for consistency with the WMAP data which finds z eq = 3141 ± 157 Q. 

For 1 + z eq = 3142(1 ± 0.05), 

7 = (1 - f3)[l - 0.503(1 ± 0.114)(3142) 3,3/2 ]. (17) 

The 68% and 95% constraints when Ob — for 7 versus j3 derived from z eq are shown in Figure [21 For (3 = 0, the 
WMAP measurement of the ear ly ISW effect Q requires 7 = 0.50 ± 0.06. 

It is easy to understand this result by comparing the high redshift evolution of the {(3, 7} model (for (3 = 0) with 
that of the ACDM model. Consistency between them requires that (1 — 7 ) 2 ~ Om = 0.26, so that 7 « 0.49, in 
excellent agreement with the more direct result from the early ISW effect. 

For each of the independent constraints from the SNIa data and the early ISW effect, the best fit occurs at (3 = 
and, as may be seen from Figure O the parameters identified by these two constraints diverge from each other as (3 
increases. In Figure [3] are shown the probability distributions for 7 when (3 = derived from the SNIa data [H and 
from the early ISW effect jij. 

The early ISW combination of parameters appears to be in some conflict with those identified by the SNIa data [H . 
For example, for (3 = and 7g NIa = 0.63 ± 0.04, leads to a prediction of z eq , (1 + z e( j)sNia = 1661^524, which is some 
2.7cr away from the WMAP value We return to a more careful discussion of the tension between the SNIa and 
WMAP data for this class of models in the context of our discussion of the more realistic model including baryons. 
Nonetheless, the results presented here provide a useful background for the subsequent discussion of those cases for 
which Ob 7^ 0. 



III. INCLUDING BARYONS 

The analysis of LSS @ established that as an alternative to dark energy, a model with the creation of cold dark 
matter by the gravitational field is capable of accounting for the evolution of a Universe in which the early-time 
decelerated expansion is succeeded by a late-time accelerating phase, consistent with the SNIa data. This encourages 
us to explore a more realistic version of the LSS model including baryons (and radiation) . While qualitatively different 
from the models explored by LSS, models with baryons are quantitatively very similar to them since for all intermediate 
and low redshifts pe <C Pcdm- As in LSS, the SNIa data Q will be used to constrain the late-time evolution of this 
class of models. To further test them and to constrain their parameters, the observational early ISW effect constraint 
on the redshift of equal matter and radiation densities inferred from WMAP, 1 + z eq = 3142 ± 157 = 3142(1 ± 0.05) [| 
is used. 

For a successful model be consistent with the early ISW effect requires that the early, high redshift evolution of 
these models be nearly identical with that of the ACDM concordance model. However, this does not guarantee that 
the growth of structure in these models need be the same as in the concordance model. For the very early evolution 
of both models only baryons (for BBN) and radiation are important; Pcdm and p\ may be neglected. In particular, 
since the results of BBN depend only on the baryon and radiation densities and are independent of pcdm and p\, as 
well as of any particle creation, it is best to choose Ob from BBN and not, for example, from observations of large 
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scale structure such as those provided by galaxy correlation functions. In the standard model, with three flavors of 
light neutrinos, the predicted primordial abundances depend only on the baryon to photon ratio which is directly 
related to Q^h 2 . To constrain f^e, deuterium is the baryometer of choice (see [Toj for a recent review and further 
references). For BBNp|, using the latest observations of deuterium [lj|, n B h 2 = 0.0218 ± 0.0011. For the HST Key 
Project value of Hq Q, = 0.042 ± 0.010. Since the models we consider are flat and contain no dark energy or a 
cosmological constant, Qcdm = 1 — Ob = 0.958 ± 0.010. These are the values adopted in our analysis below. 

A. n B + 0,/3 / 0,7 = 

As noted earlier and demonstrated in LSS, in the absence of baryons and when 7 = the expansion of the Universe 
always decelerates for < (3 < 1/3 and always accelerates for 1/3 < j3 < 1. However, with the inclusion of baryons, 
which dominate over the CDM at high redshifts, the early expansion of the Universe always decelerates. In this case, 
for a choice of (3 in the range {1/3, 1}, the late time expansion will accelerate. To better understand this transition, 
consider the "total equation of state" , w = ptot/ptot, where ptot = Pb + Pcbm + Pc = Pc, since pb = Pcdm = 0. 
For w > —1/3, the expansion decelerates, while for w < —1/3, it accelerates. Including baryons, but ignoring the 
early-time contribution from radiation, 



w = Pc/Ipb + PCDm] 



Since, for 7 = 0, p c = -/3pcDM, 

w = -(3 



1 



(—) 

\ PCDM / 



-1 -1 



= -/3 



fir 



1 - n T 



(1 + *) 



3/3 



At present, for z = 0, w = — /3(1 — J7b), so that for J1b = 0.042, the expansion is accelerating if j3 > 1/3(1 — 11b) 
In the future (a — > 00, z — > — 1), w 



(18) 
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FIG. 4: The equation of state parameter w versus redshift relations for Ob = 0.042 and 7 = for several choices of j3. Also 
shown for comparison are the w(z) variations for the Einstein-deSitter model (w = 0) and for the ACDM model (solid curves). 
Notice that z < corresponds to the future (a > ao) evolution of the Universe; as z — » —1, a — > 00. 



As may be seen from eq. 
occurs at redshift z t where, 



19, for 7 = and 7^ the transition from deceleration to acceleration 



-1/3) 



Zt 



1 - fin 
fin 



(3/3 - 1) 



1/3/3 



(20) 
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The z t — (3 relation for 7 = and VIb = 0.042 is shown in Figure [3 Without any fine-tuning the expansion of the 
Universe has a late-time, low-redshift (0 J$ z t < 4) transition from decelerating to accelerating. The present epoch is 
"selected" in this model by the BBN-determined value of f2 B 7^ 0. 




FIG. 5: The transition redshift, zt, as a function of f3 for 7 = and fie = 0.042. 



In this case (7 = 0), as may be seen from Figure [71 the best (albeit not very good) fit (see Fig. [T]) to the SNIa 
data H is for P = 0.61, so that z t = 4.0. For this model the baryons dominate over the CDM during the early 
evolution of the Universe, when z > z*, where pb(z^) = pcdm(^), 



1 + 



1 - rii 



1/3/3 



(21) 



For the model with {3 — 0.61 and f^B = 0.042 (and 7 = 0), z* = 4.5. At high redshifts, z 3> 5, pm — ► Pb, so that 
1 + z eq — > ^Ib/^Ir = 522, clearly in conflict with the constraint on z eq from the early ISW effect. So, although it is 
of some interest that the late time evolution of this model with 7 = and f^B 7^ is capable of a transition from 
decelerated to accelerated expansion without any fine-tuning, the early evolution of this model is in conflict with the 
WMAP data. 

B. n B 7^0,7 /0,/? = 

Motivated by the LSS result for J7b = that the best fit to the SNIa data occurs for (3 = 0, we next explore this 
case with J7b = 0.042. Setting (3 = but allowing for a non-zero value of the baryon density, the Friedman equation, 
including radiation, becomes 



H 
Ho 



n R (i + z) 



(1 + zf{n B + (1 - ft B )exp[3 7 (T - t )}}. 



(22) 



For the late-time evolution of H = H{z) needed for comparison with the SNIa data, the contribution from the 
radiation term may be neglected. Since VLb < 1- ^b (and, it can be checked that Q,b -C (1 — fiB)c~ 37ro ), it is a good 
approximation in solving for the z versus t relation to neglect the f2s term in eq. 19. If so, then at late-times when z 
is not too large (and r is not <C To), 



H 



(l-fiB)(l + ^) 3 exp[37(T-To)]. 



(23) 
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In this approximation, accounting for the (1 — Qb) prefactor, eq. 21 is identical in form to eq. 7 for Qb = 0, so that 
the redshift - age relation is well approximated by 

a(t)/a Q = (1 + z)- 1 « (1 - fi B ) 1/3 7" 2/3 e- 7T0 (e 3 ^ /2 - 1) 2/3 . (24) 

Evaluating this expression at t — f ( z — 0) establishes, in this approximation, the connection among t , 7, and 
(1-^b), 

cxp(-3 7 T /2) i»l-7(l- B )" 1/2 =l-i, (25) 
where 7' = 7(1 — Qb)^ 1 / 2 - The present age of the Universe (t = to, r = To) is 

r = H t = Ai n (-L-) . (26) 



37 \ 1 — 7 ' 

Solving for the age - redshift relation, 



exp(3 7 r/2) « 1 + (jz^) ^ + ^ ' ( 27 ) 

Using this approximate r vs. z relation and neglecting for the moment the contribution from the radiation density, 
the Friedman equation for H = H(z) reduces to 

W Qb(1 + Z)3 + (1 ~ r!e)(1 + z ) 3c37(T_T0) ~ "b(1 + zf + (1 - ObW + (1 - 7 ')(1 + zf' 2 ] 2 . (28) 

Notice that for z = 0, H = i/o and, for J1b = 0, this result agrees with that in eq. 12 for (3 = 0. In the high redshift 
limit, when r < ro and, including radiation, 

(|-) 2 - fi R (l + zf + [Ob + (1 - B )(1 - 7 ') 2 ](1 + zf. (29) 

If the high redshift evolution of this model is to be consistent with that of the ACDM concordance model, £Ib + (1 — 
^b)(1 - 7') 2 = 0-042 + 0.958(1 - 7 ') 2 « ^m(ACDM) = 0.26. This suggest that 7' « 0.52 and 7 « 0.51. If, instead, it 
is required that the redshift of equal matter and radiation densities agree with the WMAP determined value [H , 

Q R (i + z eq )=n B + (i-n B )(i-^') 2 , (30) 

we find 7 = 0.52 ± 0.06. For (3 = and 7 = 0.52, To = 0.97 and the present age of the Universe in this model is 
<o = 13.2 Gyr. 

As a check on our approximation to the full Friedman equation, where pb was neglected compared to pcdm in order 
to find the z — t relation, we note that as z decreases, the ratio of pb to Pcdm decreases from Ob/(1 — fiB)e _37r ° ~ 0.20 
at high z, to Ob/(1 — Ob) ~ 0-04 at z = 0. From early times to the present, the number of cold dark matter particles 
in a comoving volume increases by a factor of ~ 4.5. Notice that for (3 = 0, the early ISW effect constraint on 7 
differs little without or with baryons, 0.50 ± 0.06 versus 0.52 ± 0.06, for £!b = and £Ib = 0.042 respectively. 

A discussion of the SNIa constraint in this case ((3 = 0) is included in the exploration of the more general case 
((3^0) considered next. 

C. fl B ± 0,7/0,/? 

Based on our analyses of the fis = case which allowed for non-zero values of both (3 and 7, along with our analysis 
of the case where S^b = 0.042 and either (3 or 7 = 0, the redshift dependence of the full Friedman equation can be 
approximated as 

#■ V « «r(i + zf + n B (i + zf + [y + (i _ y _ ^(i + z fWlA 2 . (31) 



H J ' (1-/3) 

For z = 0, iT = i?Q) and the age of the Universe is 



r = H t = |-ln f - 1 / a ) . (32) 



3 7 - (3 
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FIG. 6: The age of the Universe in units of the Hubble age, to = -ffo^o, as a function of 7 for selected values of /3, for models 
with and without baryons. 




FIG. 7: As Figure [2] now including baryons. The 68% (red) and 95% (green) contours in the f3 — 7 plane, derived for Qb = 
0.042 from the SNIa data 0], are shown along with the 95% confidence band (black; see the inset which zooms in on this band) 
consistent with the early ISW constraint that z eq = 3141 ± 157 [ij. The SNIa data are best fit at (3 = and 7 = 0.66. The 
solid black curve (see eq. 34), shown in more detail in the inset, is for the best fit relation between /3 and 7 consistent with the 
observed value of z eq . 



In Figure [6] are shown the To — 7 relations for selected values of Ob and /3. For 7 > 0.45 and (3 > 0, To ^ 0.9. For 
7 = 0, t = 2/3(1 — Ob)~ 1 ^ 2 (1 — Notice that as 7 — > and (3 — > 1, t — > 00, as it should since this case is 

equivalent to the steady-state model where there is no beginning to the Universe. 



11 



At high redshifts, 

{§0) 2 * ^ r(i + z)i +nB{i+ zf " b) C~i-/ ) 2(i + zf[1 ^ ] - (33) 

Equating p^ to pm = Pb + Pcbu at z = z eq leads to a constraint on the 7 — (3 relation, 

7 = (1 - P)[(l - «b) 1/2 - (O r (1 + z eg ) - B ) 1/2 (1 + «e,) a/J/a ]. (34) 




FIG. 8: The likelihood functions for 7 inferred from the SNIa data (solid/gray) and from the z eq constraint (dot-dashed/blue) 
for P = and Q B = 0.042. 



In Figure [7] are shown the 68% and 95% contours in the (3 — 7 plane for f^B = 0.042 from the SNIa data [|[ and 
from the early ISW effect (eq. 34) constrained by the WMAP value of z eq [4] . For the SNIa data [|| the best fit occurs 
for (3 = and 7 = 0.66 ± 0.04, in some conflict with the best fit to z eq for (3 = 0, for which 7 = 0.52 ± 0.06. As may 
be seen from Figure [7J these two constraints diverge further from each other as (3 increases from zero. In Figure [5J 
are shown the probability distribution functions for 7 derived from the SNIa data [|| and from z eq Q for j3 = and 
fis = 0.042. Note that there is some tension between these two independent constraints. Indeed, if the best fit SNIa 
values of (3 and 7 are used to predict the redshift of equal matter and radiation densities, (1 + z e(J )sNia = 1798^552^ 
which is some 2.5<r away from the WMAP value [4(. 



IV. SUMMARY AND CONCLUSIONS 



As an alternative to the standard, ACDM model, we have explored here a class of models whose late-time acceleration 
is driven by the creation of cold dark matter. In constrast to the ACDM model, the models considered here have 
no cosmological constant or vacuum energy. While the dark energy models have two or more adjustable parameters 
(the dark energy equation of state, w, and its evolution with redshift, in addition to the choice of IJde), the flat 
ACDM model has only one free parameter, J7a = 1 — Ocdm- Here we have considered a class of two-parameter 
({/3, 7}) models. We found that the supernovae data prefer (3 = and that the early ISW effect is consistent with 
this choice. For this subset of creation-driven models, with only one free parameter (7), the high-redshift evolution of 
the Universe is qualitatively indistinguishable from that of the ACDM model, while the recent evolution is sufficiently 
similar to it to allow consistency with the SNIa data. From the high redshift constraint on the WMAP determined 
value of z eq Q provided by the early ISW effect, we determined 7(ISW) = 0.52. For (3 — and this value for 7, 
H to = 0.97, consistent with the estimate H t = 1.00 from the ACDM model. This choice of 7 corresponds to the 
WMAP determined value of z eq = 3142, consistent, within the uncertainties, with the ACDM value of z eq — 3223. 
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However, for ft = 0, this choice of 7 provides a poor fit to the SNIa data @. In contrast, the SNIa data prefer 7 = 0.66 
which, for j3 — 0, corresponds to H t = 1.13, in good agreement with the ACDM value and with estimates of the 
age of the Universe. But, for (3 = and 7 = 0.66, the redshift of equal matter and radiation densities is z eq = 1798, 
in conflict (at ~ 2.6<r) with the WMAP determined value. 

While these models, which are consistent with the SNIa data, offer an intriguing alternative to the standard, ACDM 
concordance model, they face challenges. For example, although the one-parameter (/3 ^ 0, 7 = 0) model with baryons 
provides a natural solution to the observed, late-time acceleration of the Universe, without any need for fine-tuning, 
its early evolution is baryon-dominated and inconsistent with the early ISW effect as constrained by the CMB data. 
In general, there is a clear tension in this class of models between the SNIa data and the independent, high redshift 
constraint from the observed early ISW effect. It should also be noted that in these models the ratio of dark matter 
(CDM plus baryons) to baryons increases during the recent evolution of the Universe from pm/pb ~ 6.0 at high 
redshifts (z > 10) to pm/pb ~ 24 at present (z = 0). Correspondingly, the baryon fraction, /b = Pb/pm decreases 
from /b ~ 0.17 at high redshifts to /b ~ 0.04 at present. While this latter value appears to be in conflict with the 
x-ray cluster baryon fraction [13| . if clusters were formed at sufficiently high redshifts, their baryon fraction may be 
representative of the earlier, higher value of /b- Due to the recent increase in the density of CDM, the late-time 
growth of structure and of the cluster baryon fraction in these models will differ from that of the ACDM concordance 
model. Therefore, before ruling out these models, it might be worthwhile to explore their late-time evolution more 
carefully, especially with regard to predictions for the CMB and for the growth of large scale structure. 
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